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Absuasi 

In this paper the effect of a small amplitude progressive wave on 
the laminar boundary layer on a semi - i nf i ni t e flat plate, due to a 
uniform supersonic freestream flow, is considered. 

The perturbation to the flow divides into two streamwise zones. 

In the first, relatively close to the leading edge of the plate, on a 
transverse scale comparable to the boundary layer thickness, the 
perturbation flow is described by a form of the unsteady linearised 
compressible boundary layer equations. In the freestream, this component 
of flow is governed by the wave equation, the solution of which provides 
the outer velocity conditions for the boundary layer. This system is 
solved numerically, and also the asymptotic structure in the far 
downstream limit is studied. This reveals a breakdown and a subsequent 
second streamwise zone, where the flow disturbance is predominantly 
inviscid. The two zones are shown to match in a proper asymptotic sense. 


*Work funded under Space Act Agreement C99066G; presently at 
Department of Mathematics, University of Manchester, Manchester, 
England, 
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1 . Int roduct ion 


Much work has been carried out on the response of a t wo -d i mens i ona I 
imcompressible laminar boundary layer on a semi -inf ini te flat plate to 
time-harmonic oscillatory perturbations (particularly those of infinite 
wavelength in the streamwise direction) of the f reest ream . Host of 

the effort has been involved with small amplitude perturbations; included 
in this category are th* works of Lighthill (1954), Lam and Rott (1960) and 
Ackerberg and Phillips (1972). Close to the leading edge, the flow is of 
Blasius type, whilst far downstream the boundary layer takes on a double 
structure, comprising an inner Stokes-type layer, and an outer Blasius-type 
1 aye r . 

Numerical investigations to this small amplitude problem, describing 
the flow from the leading edge, to far downstream have been conducted 
by Ackerberg and Phillips (1970) and Goldstein et al (1983). Some details 
of the far downstream behaviour of this problem have in the past been 
the subject of some controversy. At large distances downstream of the 
leading edge, a set of e i genso 1 ut i ons must exist, which decay exponentially 
fast downstream in the streamwise direction (this reflects the effect of 
the particular conditions prevailing upstream). One set of e i genso 1 ut i ons 
found by Lam and Rott (1960) and by Ackerberg and Phillips (1972) is 
determined by conditions close to the wall, and has decay rates that 
decrease as the eigenvalue increases. A second set of e i genso 1 ut i ons 
has been found by Brown and Stewartson (1973a, b) , and these are determined 
primarily by conditions at the outer edge of the boundary layer, and 
are characterised by a decay rate that increases with increasing order. 

Some discussion of this apparent discrepency is to be found in 
Goldstein et al. (1983). 

Further downstream still, Goldstein (1983) showed how the e i genso 1 ut i ons 
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of Lam and Rott (1960), which initially decay, develop into the 
classical, high Reynolds number limit for To 1 1 mi en-Sch 1 i cht i ng modes 
(including the unstable mode) for Blasius-type boundary layers. 

Although the Lam and Rott (1960) eigensolut ions decay downstream, they 
also oscillate with increasingly rapid (spatial) frequency, and far 
downstream, transverse and streamwise gradients must become 
comparable, giving rise to an entirely new structure, in which transverse 
pressure gradients play a key role. 

Returning to the boundary layer problem, non-small oscillations 
of the freestream (but still limited to non-reversing freestreams) have 
been tackled by Moore (1951), (1957) and Pedley (1972) for regions close 
to the leading edge of the plate, and it is again found that the flow is 
Blasius like. Lin (1956), Gibson (1957) and Pedley (1972) studied the 
far downstream region, and found the flow takes on a double structure 
similar to that observed in the analogous small amplitude case. 

Duck (1989) has presented numerical results which extend from the leading 
edge, to far downstream, in this case. 

All these aforementioned papers were concerned with purely temporal 
flow oscillations (i.e. oscillations of infinite streamwise wavelength). 
The effect of a small amplitude progressive wave on an incompressible 
boundary layer has been investigated by Rest in, Maeder and Wang (1961) 
and Patel (1975). The former authors considered the low frequency limit 
to the problem, for the very particular case when the wavespeed equalled 
the mean far field velocity. Patel (1975) gave results for both high 
and low frequency, the latter being obtained using an empirical momentum 
integral appraoch. A number of experimental results were also presented, 
and compared with theory. 

In this paper we consider the effect of a two-dimensional small 
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amplitude progressive wave in a supersonic freestream, on a laminar 
boundary layer on a flat plate. Here the situation is rather different 
to the incompressible case in a number of respects. Although eigenfunctions 
similar to those of Lam and Rott (1960) undoubtedly exist, these are not 
expected to develop into growing Tol lmi en-Schl i cht ing waves (at least 
not in the case of two-dimensional flows - see Ryzhov and Zhuk 1980 and 
Duck 1985). Although planar To 1 1 mi en-Schl i cht i ng waves are present 
downstream, these all decay and so are of little consequence. In this 
paper, as well as presenting results for the boundary layer region, we go 
on to consider the development of the flow far downstream, which turns out 
to be predominantly inviscid. A further difference found is that in the 
case of these supersonic flows (unlike incompressible flows) if the 
progressive wave is to satisfy the appropriate governing equations in 
the farfield, the spatial and temporal wavelengths are linked. (A somewhat 
related study, of a forced flow but at finite Reynolds numbers has been 
made by Mack 1975) . 

The layout of the paper is as follows: in Section 2 we consider the 
"fore region*' in which the flow close to the plate is governed by a form 
of the compressible boundary layer equations. We look at the downstream 
limit of this zone, and then in Section 3 we consider the downstream 
region, wherein transverse variations of pressure become important; 
the regions studied separately in Sections 2 and 3 are seen to match 
(in an asymptotic sense). In Section 4 we show how results from the 
previous two sections correctly match, and finally in section 5 we 
present a few conclusions. 
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2. The fore region 


We consider the viscous supersonic two-dimensional flow over a 
semi -infinite, stationary plate. We take the origin at the leading edge of 
the plate, x* the coordinate along the plate, y* the coordinate normal 
to the plate, and u* and v* the velocity components in the x* and y* 
directions respectively, t* denotes time, p* the fluid density p* 
the viscosity of the fluid, p* and T* are defined to be pressure 
and temperature in the fluid; Cp(c v ) is the specific heat of the fluid 
at constant pressure (volume), and a - M*Cp/K* is the Prandt 1 number, 

(K* being the coefficient of thermal di f fus i vi ty) ; Cp, c v and <j are 
all assumed to be constants. Subcript wi 1 1 be used to denote 

farfield, unperturbed conditions. 

The fluid is also assumed to satisfy the following equation of state 

p* - p* R*T* (2.1) 

where R* - Cp - c v . (2.2) 

We consider here, first the fluid in the far field, where the effects 
of viscosity are expected to be negligible. In this region, we take the 
flow to comprise a uniform steady stream parallel to the plate, perturbed 
by a small amplitude (0(c) « 1) progressive wave. 

We therefore write 

u* - U* + e u*(x*,y*,t*)+. . . 

v* - £v*(x*,y*, t*)+. . . 

P* p* + epi*(x*,y*,t*)+... 

T* - T* + cT*(x*,y*,t*)+. . . 
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The leading order equation of state gives 


★ * ★* 
Pcc * Pcc RToc ■ 


(2.4) 


Subst i t ut ion of 
equation yields (at 


(2.3) into the momentum equations and continuity 

0(e)) 


r _ * ^ 

* OUT , .** OUT 

p “ 87 * + u “ 



(2.5) 


* 

Pcc 



+ 


u, 


dx* 



( 2 . 6 ) 




+ 



* 9uT 

+ Poe r* 

dy 


0. 


(2.7) 


Differentiating (2.5) with respect to x*. (2.6) with respect to y* , 

and adding yields 




& + + ' 


( 2 . 8 ) 


If we write 

a x 2 - 7R*T tt * (2.9) 

with > - c p/ c v (assumed constant) (2.10) 

where a a denotes the speed of sound in the far field, then (2.8) 

* 

may be written entirely in terms of pj , viz 


*2 



( 2 . 11 ) 


We now seek progressive wave solutions of this equation, writing 

p* - p*jexp[ia*(x*-c*t*+xV)] + exp[ io*(x*-c*t*-X*y*) ]} (2.12) 
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(corresponding to an oblique wave, in general), where p^* is a constant. 
Substitution of (2.12) into (2.11) reveals that we must have 


c 


* 


* 




1 i ^ (l+x* 2 ) i . 


(2.13) 


where denotes the unperturbed freestream Mach number, 

It is now possible to write down the progressive wave solution for 

* * 

U} and vj > namely 


* 

U 1 - 


*jexp [ ict*(x*-c*t*+X*y*) ] + exp[ ia*(x*-c*t*-X*y*) ] 

(U«-C*)Pecl 


(2.14a) 


* Xpi 

v ! “ - - ±± 


. ,exp( ict*(x*-c*t*+X*y*) - exp [ iQ*(x*-c ir t ir -X*y 1,: ) ] } . (2 . 14b) 


'k / * 


Equation (2.14b) also satisfies the normal flow condition, although 
the no-slip condition is violated by (2.14a), since 


ui(y*-0) 


1 eL 


exp{ iQ*(x*-c*t*)] . (2.15) 


(U* _c *)P<x 


Notice in particular that in the case of plane waves (X* - 0) , 


u oc [i i 


(2.16) 


* ^ 

vi - 0, 


(2.17) 


(whilst (2.15) is unaltered). 

The slip velocity is reduced to zero, in the usual way, by the 
inclusion of a thin boundary layer. The boundary layer approximation 
reduces the governing equations to (Stewartson 1964) 


3u* * Bu* * 3u* 

— T + u — - + v — x 

9t* 8x* By* 


_ . V + _L 


Bx* ' By* ^ By* ' ’ 


( 2 . 18 ) 


7 


1$ * S' ( '* u * ) * £ < '* v * > - °- <2 ,9) 


* 

3t* + * 3T* . * 3T*1 

dp* * 3p* a 

* C P 3T *‘ 

* 

3u* 

P c p 

[at* 3x* 3y*j 


* ° 3y* 

+ P 

|ay*J 


( 2 . 20 ) 


1 ^- 0 , ( 2 . 21 ) 

3y 

(together with (2.1)). 

These equations can be simplified by the use of a generalisation 
of the Howart h-Dorodni tsyn transformation, suggested by Stewartson (1951) 
and Moore (1951). We write 


P 


* 


* 9y 

Poo 


ay* 


* 


U 


3y 


( 2 . 22 ) 


where Y now replaces y* as the independent transverse variable, and 
v* is essentially a streamfunct ion. Then 


* 

Ez. 

p* 


lav*' 

+ u* 

a? ' 

+ 

a? ' 

l 

.[ax* 

Y,t 

3x* 

y, t 

at* 



(2.23) 


and the equations of motion and energy reduce to 

aV + B4* ay _ av 
a?at* a? ax*a? a? 2 ax* 


a_ d 2 v*l 

a? a? 2 J ’ 


T 

* * 
T 1P1 


§E1 

3x* 


(2.24) 
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a r 

8t * 


* * 
PcoPI 
* * 
PooPoo 


8t* 

By:* 8T* 

T* 

[dp* . d v '* 

- *1 

opi 

8x* 

8x* 8Y 

* * 
C P P1 T 1 

iat* 

8y 

8\*. 

a_ 

U*t » 

★ JL, * 

, PIP T® 

8 2 ^*' 


a? 

o/^t* a? . 

* * * 
c pPooP®^ 

a?2 . 



(2.25) 


On account of (2.21), we have written 

p*(x*,Y,t*) - p^x*,!*). (2.26) 

The usage of the streamfunct ion ensures that the continuity condition 
(2.19) is satisfied. 

We now seek a solution to (2.24) - (2.25), subject to the boundary 
condi t ions 

— (Y - 0) - 0, (2.27) 

a? 

T*(Y -»«)-» 1 + 0(e), (2.28) 

^*(Y - 0) - (Y - 0) - 0, (2.29) 

Y 

* * 

\f* -» U* - — * exp[ io*(x*-c*t*) ] as Y (2.30) 

Y (Uee-C *) Peo 


At this stage we must specify the particular viscosity/temperature 
model, and here we choose the simplest example, namely the linear law of 
Chapman (See Stewartson 1964), 




* 



(2.31) 


We now seek a perturbation solution (in powers of e) ; we write 
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4* - + c ^* + - 

T* - Tq + «T* +. 

p* - P0 + ep*+ • • • 


Y - Y 0 + cY + . . 

* ★* — * 

Pi - PooRToo + CP +. 


(2.32) 


Substituting (2.32) into (2 . 24) -(2 . 25) yields, to 0(c®) 


3^0 3 2 \tp _ 3i/-p 3 2 ^q _ p® 3 3 v'q 

3Yq 3x*0Yq 3x* 3Yq p® 3Yq 


(2.33) 


3^n 3 Tq _ 3^o 9lfl _ j£_ 9 2t Q + 3 ^q 

dYp 3 x* dx* <5 Yq * iy2 * -ly^ 

u u ffp® 3 Yq e-pa, IoYq 


(2.34) 


* ★ 

The solution of 4o (and Tq) is now routine, and merely 
corresponds to the similarity form of Blasius. Writing 


4o ~ 


*-* .*1 i 


2pcoU®x 


* 

Pec 


F 0 (iJ>. 


(2.35) 


Tq “ T ® G 0 (tj) , 


(2.36) 


where rj - Yq 


* * ii 

Pcr^cc 


2p®x’ 


(2.37) 


then 


F 0 + F 0 F 0 - 0- 


(2.38) 


\ Gq + F 0 C 0 + (7”1)Mo 5 2f 0 “ 


(2.39) 
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The boundary conditions to be applied are 


F 0 (0) - Fq(0) - 0, 

F 0 (») " 1, 

Gq(0) - 0. 

G 0 (®) - 1 . 

The 0(c) equations are rather more complicated, although 

derivation is routine. The momentum and energy equations at tl 
are, respectively 



**Y 0 





“Ox 


-7* y * r* 

* ^Y 0 V 0 ‘ V X* Vo Y 0 


Tn 9p* , £ a 3 ?* , P*^I e 3 4 

TooPoc 9x* Poc 9 Y 0 3 PcoPoc 3Yq 3 


* 

T 

t* 







I ★ 

v x* 


* 


. T o 
* * 
CpPocTcc 


Bp! 

at* 


+ 



a? !' 

8x* 


* 


* 

Pccff 




* 

PwO 


T 


oW 


(2.40) 

their 
is order 


(2.41) 


(2.42) 
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Our primary concern here will be with the momentum equation (2.41) 
(although a similar treatment may be carried out on (2.42)). 

We already have (2.35)— (2.37) t and now we write 


iPl* 


★ 4. * 

(Uoc-C*)P« 


2/^x* 
* * 
Poo ^ 


i 


Fi(t hi) 


-ia*c*t* 


(2.43) 


where £ - x*$. (2.44) 

These transformations, when applied to (2.41) yield the following 
equation for Fj(ij,£) f 


2 i Ot c* 

+ Fq^I^^ + ^l^Orjrj + ^ ^Or? ^1 tj£ + ^1£ ^Orjrj 

tyr 


2Gq(t?) £ 2 ia 

[-H 

. 1£2 “* - -A. 2 

[-£->] 

i£ 2 a* 

e ^OrjrjTj * 

(2.45) 


lu c J 


lu a J 




a* may be completely eliminated from the problem, by introducing the 
new scaled variable £ 


l 


★ ★ 

c or 


(2.46) 


It is also convenient to introduce the non-dimensional wavespeed, 

c - c */u£ . (2.47) 

and so (2.41) becomes 

+ + F 0tj»j + ” ^Ori f 1tjj + 

“ 2 iC o (t?) S 2 (1- \) e^ 2/c - *^2 (c-1) e 1 * 2 ^ F 0lJ1J)J , (2.48) 
subject to 

Fi(r~0) - Fi^(ij-O) - 0 


as t; -» <* . 
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(2.49) 


This equation may be viewed as the compressible counterpart of the 
work of Patel (1975), and also of Ackerberg and Phillips (1972) (and others 
cited previously) which was specifically concerned with temporal 
osci 1 lat ions. 

Equations (2 . 38) -(2 . 39) were solved numerically, using a second-order 
finite-difference scheme, in which the momentum and energy equations were 
split into three and two first order ordinary differential equations 
respectively, and solved iteratively using Newton's method. Fq(t;) and 
Cq(t?) were then input into (2.48); this (linear) equation was solved using 
second order finite-differences in 17 , Crank-Nicolson marching in £, 
again by splitting the equation into a system of first order equations 
(in 7 ?). 

— i 1 2 /q ^ 

Results for e F 1^(77 - 0)/£ are shown in Figs 1-4. In all 

cases we set 7 * 1.4, o — 0.72, X - 0 (plane waves) and the results 
are believed to be correct to within the graphical tolerance of the 

figures. Figures la,b show results for * jl , c - 1 - Figs. 2 a,b 

Moo 

for M* - Jl, c - 1 + i/Hoo, Figs 3 a,b for - 5. c - 1 - l / Hx , 

Fig. 4a, b for M* - 5 , c - 1 + Vmcc- 

It is readily apparent from these results that as $ -> » (i.e. far 

-i£ 2 / 

downstream), e /c F^ ^(rp-OVf asymptotes to a constant, and it is 
a simple matter to confirm this analytically. 

From (2.48), we have as £ o° 

2il 2 F lrj - hovFlnl + o,, - 2 iGqCtj) (1 _ l) e ,£ 2 / c. (2.50) 

The inclusion of the second and third terms here is most easily 
justified by transforming to the (non-dimensional) streamwise (x) 
coordinate, where 

x - £2 f (2.51) 
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whence 


*^1 rj ~ ^Otj ^1t;x + F 1x " JCq(i?)( 1 - “)e . 

Writing 

Fj(t?,x) - F 1 (»})e iX//c , 

then 

Fld> - (C-F 0 ,) |” ' 

0 

This solution generally violates the no-slip condition 

Fi^rr 0 ) - C o( 0) (1-1). 

However this is easily rectified by the inclusion of a 
where 


y - - V xi - 0 ( 1 ), 

and 

F 1Y yy + 2iF 1Y - 2iC 0 (0) (1 - l)e ix/c . 

Writing 

Fj(Y,X) - fi<Y) e ix/c > 

then 

flY <Y) - G 0 (0) (1 - 1) |l - e' (1 " OY j 

which does satisfy the no-slip condition on Y-0. Equation 
then gives 

F ln ,OrO) - ie if2/c G 0 ( 0 ) (1 - 1) (1-i). 

This is shown as an asymptote on Figs. 1-4, and the result 
confirmed. Equation (2.54) does admit the possibility of a 


(2.52) 

(2.53) 

(2.54) 

on t^O, since 

(2.55) 

Stokes layer, 

(2.56) 

(2.57) 

(2.58) 

(2.59) 
(2.59) 

(2.60) 

is seen to be 
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singularity, wherever 


Fo V ~ 1 / c . (2.61) 

and a critical layer is necessary if this occurs, as described in the 
Appendix. 


Unfortunately the above expansion is not uniformly valid as £ -> oo. 

This is most easily seen by comparing the magnitude of the transverse 

► 

pressure gradient (3p* . ) with the transverse inertia terms (in particular 

/ ay 

the p*u*v* term) namely 


V 

ay* 


o 


* *2 
Poo Uoo 


* * 
P m U„ 

oo w oo 


P CO X 



* * 4.1 

Moo X 


P® <** 2 U^ 

★ 

P® 



- 0 


Poo ^00 

M^x*a * 2 


(2.62) 


Hence there will be a breakdown to the above solution when 


x 


* 


0 


* * 
Poo U® 


★ *2 
Moo <* 


(2.63) 


when the transverse pressure gradient can no longer be neglected. In the 
following section we go on to consider the repercusions of this. 

One final point is that eigensolut ions analogous to those described 
by Lam and Rott (1960) are certainly expected to occur far downstream 
(although their magnitude was such that they were undetectable in our 
numerical results); however, as noted in the previous section, unlike the 
subsonic case, in the supersonic planar case currently under cons i derat i on , 
these eigensolut ions are of little consequence. 
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3. The far downstream region 


As noted at the end of the last section, unsteady transverse pressure 
fluctuations will become significant in the physics of the boundary layer 
when 


★ 

* 

Pod 

Uoo 

* 

*2 


a 


This implies (and is consistent with) 

y* - 0(a*-l), 


(3.1) 


(3.2) 


implying the transverse scale (i.e. the boundary layer thickness) is 
comparable with the wavelength of the imposed wave. At the same time, the 
streamw ise lengthscale of the wave (a*~l) is much shorter than the 
developmental length of the boundary layer (i.e. (3.1)), and so the 
parallel flow approximation is a rational procedure, and is adopted 
here . 

We choose to make the following non-dimensional isat ions 

y - 

a — a* 6* , 


(U,v) - (u*,v*)/U*, 

* * 

P “ P /Poo. 

. it ^ it 

P - P*/(PooR*Tco) . 

T - T*/Tco 

C - c*/U* , 

it * 

t - 

< 5 * 


( 3 . 3 ) 


Here 6* represents the boundary layer thickness, (and is the key 

lengthscale used in the non-dimensional i sat i on process), i.e. 
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We then expect the solution to develop as follows 


u - U(y) + ef (y) E+. . . 

v - ea *>(y)E+. . . 
p - 1 + e x(y)€+. . . 

T - T(y) + c 0(y) E+... 
p - p(y) + £ R(y) E+. . . 


with E - e 


ia(x-ct ) 


(3.5) 


and where U(y) , p(y) and f(y) denote the mean flow speed, density 
density and temperature respectively. 

The governing equations are then found to be (Lees and Lin 1946) 


p[ i(U-c) f + UVI - ^2 

t 

pio 2 (U-c),p - ^£2 
i(U-c)R + p(y?'+ if) + p' <p - 0, 
r e 

r - - + - , 

p i 


(3.6) 

(3.7) 

(3.8) 

(3.9) 


[ i (U-c) 6 + f' yj] - 'Prl t(U-c)T (3.10) 

After some simple manipulation (following Lees and Lin 1946), 
we arrive at the following equation for <p, 


d_ [ (u-c)g'-uy 

dy It-m£(U-c ) 2 . 


a 2 (U-c)vp _ 0 

T 


(3.11) 


(It is worth noting at this stage that the analogous governing equation 

of Mack 1984, 1987 appears inconsistent with this equation, and should be, 
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In Mack's notation 


(a U-u) Dv - a DUv 

c? (a U-co) v 

(1-M 2 )T 

f 


Finally, we follow the same approach as used previously, by using 
the trans format ion 


(analogous to 


► 


By 

3y 


p (y) 


(2.21)) and so (3.11) becomes 


(3.12) 


d_ 

dY 


[U (Y) - c] ^ - u’(Y)*> 
T[T-Moc 2 (U(Y)-c) 2 ] 


-a 2 (U-c)v? 


0 , 


(3.13) 


where U(Y) - Fq(»?), 
Y - V. 

T(Y) - C 0 (ij). 


(3.14) 


with Fq(tj) defined by (2.37), and Gq(t?) 
The boundary condition on Y - 0 is 


by (2.38). 


yp(Y - 0) - 0, 


(3.15) 


corresponding to the impermeability condition. As Y we require 

the match with the progressive wave solution. From the previous 
section, this condition can be described by the superposition of an 
incoming and outgoing wave, if 


ioXY ^ 
tp — e + 


_ -ioXY 
0e as 


(3.16) 


0 is a reflection coefficient, which is to be determined as part of 
the solution. The wavespeed c is given by 
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c 


(3.17) 


-111- (i+x 2 )i 


on account of (2.13). 

The problem then may be thought of as that to determine 
X (generally complex), Moo, and a* This was achieved by 




^(Y) + e 


iaXY 


where , 


fU(Y)-c] g (Y) - 0’ (Y)y? 

dY T(T) [T(T)-M 2 (U(T)-c) 2 ] 
00 


a2(U(Y)-c)^(Y) 


iaX Y 

- e 


q 2 X 2 (U-c) + u"(Y ) 

T[T-m£(0-c) 2 ] 


+ 11 < Y) 

T 2 


iaX(U-c) -O' (Y) 
[T-m£(U-c)2] 


+ [iaX(0-c)-u'(Y)][f’-2M^(U-c)0'(Y)] 
T[T-m£(U-c ) 2 ] 2 


This equation was split into two first order equations, 


/3, given 
writing 

(3.18) 


(3.19) 

one for 
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<?' (Y) , the other for F' (Y) , where 


F (Y) 


[ 0(Y)-cjg'(Y)-U' (Y)£ 
T(T-m£(U(Y)-c)2] 


(3.20) 


A simple fourth-order Runge-Kutta method was then used. To obtain the 
appropriate behaviour at infinity, the radiation condition 


p + iaX^- 0 (3.21) 

Y 

was applied at some suitably large value of Y. 

This solution will in general fail to satisfy the no-slip condition; 
however this is easily remedied by the presence of a Stokes layer of 
thickness Y - 0(1) (see (2.56)). 

Figures 5-10 show the variation of the reflection coefficient 0 
(see (3.12)) with a. Results for - jl , X - ^ » negative root 

of (3.17) are shown in Figs. 5a, b; for - j2 , X - ^ , 

negative root of (3.17) in Figs. 6 a, b; for - /2 X — ^ + ^i f positive 

root of (3.17) in Figs. 7a, b; for * j2 , X - ^ + ^i, positive root of 
of (3.17) in Figs. 8 a, b; for - 5 , X - j + i, negative root of (3.17) 
in Fig. 9a, b; for M^-5, X - g + g i, negative root of (3.17) in 

Figs. 10a, b. All these results were obtained with 7 - 1.4, a - 0.72. 

The primary result which is clearly visible from these computations 
is that 0 -> -1 as a 0. This indicates a correct match with the 
results of the previous section in the upstream limit (a 0 , here, 
may be interpreted as the upstream limit), in particular with (2.15). 

This match will be studied in full in the following section. 

Lastly, it is also clear that as a -» « (downstream limit), 0 -» 0 . 

In the following section we summarise the conclusions of this work. 
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4. The matching of the fore and downstream re&Lons 


Let us first consider the limit of (3.11) as a -* 0 (which as 
pointed out at the end of the previous section is equivalent to the 
upstream limit of the downstream zone). Two lengthscales for y 
emerge, the first, where 

y - ary — 0(1) (4.1) 


involves a solution to (3.11) of an incoming and outgoing 
wave , name 1 y 

v , - e iXy + 0e' 1Xy + 0(o°). (4.2) 


We may impose the impermeability constraint on this system, which 
requires 0 — -1 and so to leading order 


<P - 


. ix >' 



(4.3) 


this condition matching correctly with the outer fore-region solution 

(2.14). 

The second important lengthscale is y - 0(1) itself, where the 
variation of T and 0 must be taken into consideration. To leading 
order we have 

d_ f (U-c) - 5’ A _ Q (4.4) 

dy l T - Moo^U-c) 2 1 


The solution to this is 




C(U-c) 


y 

o 


T-M 0C (U-c) 2 

(0- c) 2 


dy. 


(4.5) 


where the impermeability condition has been imposed, and C is a 
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constant which is determined by matching with (4.3), which yields 


C 


2cri (1-c) 
X 


(4.6) 


This completes the details of the a 0 solution; however for 
the purposes of matching with the fore-region, the simplest illustration 
is by means of the perturbation streamwise velocity f(y). By (3.6) 
we see we may write 


f(v N 


-tt toy 

(U-c) (U-c) 


(4.7) 


The pressure term x is (by (3.7)) seen to be independent of y, and 
given by 

* - -iC-yM* 2 . (4.8) 


and so 


f(y) 


iCT 

U-c 


+ iU C 


y 

o 


T 

(U-c) 2 


dy 


-iU’ CM^y, 


(4.9) 


(if 0 - c, then a critical layer of the type described in the 
Appendix will be present). 

We now match the downstream limit of the fore-region w ith (4.9) 

The 0(e) perturbation to the x* component of velocity, as 
x* « may be written 

' * ic**(x*-c*t*) - * 

=* <**•"> - ~ . - F i= + "- %}<*>• < 41 °) 
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where rj Is defined by 


V ~ Y 


★ * 

Poo Ug 


2poo X x 


*v* 


C4.11) 


whilst 


- 


’ * - 

p ( ^-pQ ( Jl d ^, 

PC v) 


(4.12) 


and F^(tj) is given by (2.54) (with tj replaced by rj) . Notice 
that 


V “ 


Pcc 


dy. 


(4.13) 


Equation (4.12) may be written 


v(y) 


p* g-] (y)+fl?*(y) 


P O c*R*T oo *C 0 (y) T OD 1c Go(y) 


r 


dy. 


(4.14) 


Here *!*(y) - ^T* C 0y (y)G 0 (y) , (4.15) 

whilst ^ 2 (y) is determined from the far downstream limit of (2.42) 
which yields 


e * (y) . c Q<y>P* + T = c Qy F l<yI E *£l, 


(4.16) 


c nPoo 


IT * 

Uq -c 


where we have written 


_* ia*(x*-c*t T ) 
t - e 


(4.17) 
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u 0 *(y) - u«*u 0 (y) - u* Fq(^), 


(4.18) 


Cl — * 

(U C0 *-c*)p 0C 


(4.19) 


Consequently we see that 


v(y) 


y 

-It; 


p* 


p OB R*T eo C 0 (y) 


^(y)Cnv(y) 

c 0 


— * 

p 


^"ooCpPooCo ( y ) 


_ Cq v Fi (y) 


c 0 (u5 


-c*) J 


The solution of this equation is 


(4.20) 


T?(y) - 


Co(y) 


p* 




£l_ . C 0v F 1 (y)E*C 1 ^„ 


★ * 


Poo^ To, U 0 -c 


(4.21) 


We now see that as x -* ® 


u*(x*,y) 


- CjEFiyCQ + l >0 (y)E' 


■ {[-£ 


CnvFi(y)c 


7P«T®R Uq-c 


^<*y. 


(4.22) 


where we have used the condition 


P - 2p E . 


(4.23) 


Equation (4.22) may now be written 


=*(x*.y) - ClE ef8^(y>°~ 2 V* > f -- - CjE*^ <y)(c-l) | 

TPocTodR 


(4.24) 


- c l E * G Q ( y) (£ll> _ CjM^EVoCyH 1 ^) - C 1 E*Uo(y) (c-1 ) 

c-u 0 (y) 


c 0 (y) 

(c-U 0 ) 2 
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dy . 


(4.25) 



(4.26) 


Sett ing Ci - 
and noting that 


we see that (4.25) 


iU*C/(l-c), 


c 0 (y) - T(y) . 

(4.27) 

uo<y) - 0(y). 

(4.28) 

matches correctly with (4.9). 
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5 . Discussion and conclusions 


In this paper a description of the effect of a small amplitude 
progressive wave on a supersonic viscous flow over a semi -inf ini te flat 
plate has been given. The original motivation in this work was to study 
the possibility of a supersonic analogy of the so-called receptivity 
problem; however it turns out that in this case there is no receptivity, at 
least in the absence of boundary curvature. However, it would certainly 
be possible that a higher order analysis incorporating boundary- layer 
growth terms (i.e. non-parallel effects) would give an element of 
boundary curvature. 

Associated with the present study, the form of the compressible 
Stokes layer has been described, together with the ultimate breakdown 
of the (boundary layer) structure of the perturbation solution which far 
downstream becomes predominantly inviscid. The analysis, together with 
our numerical results, formally indicate a proper match between the two 
regimes . 
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Appendix 


Here we consider the critical layer which will occur if 

F 0 ?7 " c (A.l) 

in (2.54). Since 0 < Fq^ < 1, this situation is only possible if the 
negative root of (2.12) is chosen. 

The treatment of the critical layer is fairly standard, and follows 
closely the structure found in, for example, linear stability theory 
(see for example Rosenhead 1963); as a consequence we shall keep details 
to a mini mum . 

Suppose that we have a line, t; - tj c , such that 

F oij (»r’Jc> “ c ; < A - 2 > 

we then write 

Fl(T),x) - e iX ^ 2 Fi(jj), (A. 3) 

and then taking the inviscid equation (2.52), and performing a 
Frobenius expansion about rj - tj c , gives, for rj < r] c 

P l(T?) _ C(i? n )(c-1) + (c-1) 0n(t7 c-v) Igo(’Jc) - F 0 ( t ?c) c q( t ?c> 1 

F 0(V C ' F 0 (,) ^ F 0 (»? c ) i 

+ 0 (i] c -ij) . (A. 4) 


For tj > rj c the same procedure yields 


Fj(,) - Cd g )(c-1) + d-l>_ ( , 1 _, )|Cn „ c _, l!l „ 
Fq^c) Fq(ij c ) 


Co(’Jc) 


_ Fp" ’ (T)r)Co(T? r ) | 

Fo(’Jc) 1 


+ 0(r) c -rj) . 


(A. 5) 


The ambiguity in sign here must be resolved by incorporating the effects 
of viscosity inside the critical layer. The similarity with linear 
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incompressible stability theory is immediately apparent, the chief 

* 

difference being the inclusion of the Go(tj c ) term, i.e. the temperature 
gradient term, which is obviously absent in incompressible situations. 
Following classical stability theory, we expect that inside the critical 
layer the solution develops as follows (assuming x»l) 


Fi(tj.x) 



FlO<0 


+ x‘ 1/3 F n (f) 



(A. 6) 


l/o 

where f - (rj c -t])x 3 

- 0(1) inside the critical layer. 


(A. 7) 


The governing equation for F]^(f) is found to be (from (2.48)) 

~.ti oj it ' • 2i 

Fio " — F 0 (i} c )fF 10 + — F 0 (iJ c )F 10 


2iC 0 (ij c ) (1 - -) 


(A. 8) 


The solution to (A. 8), which matches on to both (A. 4) and (A. 5) is 
simply 


r 10 - C ° ( V (C - 1) - 

Fo(’lc) 


(A. 9) 


The equation for Fji(f) is found to be 


Fii - Fo(Tj c )rF n + f"o 03c) F 11 “ 2iC 0 (i» c )r(l - 1 -) - 2 ± F 10 fFo 0j c > • 


(A. 10) 


The solution of this may be expressed in terms of the functions due 
to Holstein (1950) (alternatively in terms of Airy functions) and this 
demands that if a proper matching is to be achieved with (A. 4) as 
f cc, and (A. 5) as f -> then we must take the negative sign in 
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(A. 5). 


Although the temperature perturbation T* exhibits an algebraic 
singularity as 17 -» ij c (with T* proportional to (ij-ijp)* 1 ), 
nonetheless viscous effects in the f - 0 ( 1 ) layer are sufficient to 
alleviate this singularity. 

This then completes the determination of the critical layer, and 
resolves the sign ambiguity in (A. 5). 

Acknowledgements 

A number of the computations were carried out at the University 
of Manchester Regional Computer Centre with computer time provided 
under S.E.R.C. Grant No. GR/E 25702. 


29 



References 


1. Ackerberg, R.C.; and Phillips, J.H.: The Unsteady Laminar 

Boundary Layer on a Semi-Infinite Flat Plate Due to Small 
Fluctuations in the Magnitude of the Free-Stream Velocity. 

J. Fluid Mech., vol. 51, pt. 1, Jan. 11, 1972, pp. 137-158. 

2. Brown, S.N.; and Stewartson, K. : On the Propagation of 

Disturbances in a Laminar Boundary Layer, I. Proc. Cambridge 
Philos. Soc., vol. 73, pt. 2, Mar. 1973, pp. 493-503. 

3. Brown, S.N.; and Stewartson, K. : On the Propagation of 

Disturbances in a Laminar Boundary Layer, II. Proc. Cambridge 
Philos. Soc., vol. 73, pt. 2, Mar. 1973, pp. 505-514. 

4. Duck, P.W.: Laminar Flow Over Unsteady Humps: The Formation of 

Waves. J. Fluid Mech., vol. 160, Nov. 1985, pp. 465-498. 

5. Duck, P.W.: A Note on a Numerical Method for Treating 

Time-Periodic Boundary Layers. Submitted to J. Fluid. Mech., 
1989. 

6. Gibson, W.E.: Unsteady Laminar Boundary Layers. Ph.D. 

Dissertation, Massachusetts Institute of Technology, 1957. 

7. Goldstein, M.E.: The Evolution of Tollmien-Schlichting Waves 

Near a Leading Edge. J. Fluid Mech., vol. 127, Feb. 1983, 
pp. 59-82. 

8. Goldstein, M.E.; Sockol, P.M.; and Sanz, J.: The Evolution of 

Tollmien-Schlichting Waves Near a Leading Edge. Part 2. 
Numerical Determination of Amplitudes. J. Fluid. Mech., vol. 
129, Apr. 1983, pp. 443-454. 

9. Holstein, H.: Uber Die Aussere und Innere Reibungsschicht bei 

Storungen Laminarer Stromungen. Z. Angew. Math. Mech., vol. 30, 
1950, pp. 25-49. 


30 



10. Kestin, J.; Maeder, P.F.; and Wang, H.E.: On Boundary Layers 

Associated With Oscillating Streams. Appl. Sci. Res., vol. A10, 
no. 1, 1961, pp. 1-22. 

11. Lam, S.H.; and Rott, N. : Theory of Linearized Time-Dependent 

Boundary Layers. (Avail. NTIS, AD-246519.) AFOSR-TN-60-1100 , 

1960. 

12. Lees, L. ; and Lin, C.C.: Investigation of the Stability of the 

Laminar Boundary Layer in a Compressible Fluid. NACA TN-1115, 
1946. 

13. Lighthill, M.J.: The Response of Laminar Skin Friction and Heat 

Transfer to Fluctuations in the Stream Velocity. Proc. Roy. Soc. 
London Series A, vol. 224, no. 1156, June 9, 1954, pp. 1-23. 

14. Lin, C.C.: Motion in the Boundary Layer with a Rapidly 

Oscillating External Flow. IX Congress International de 
Mecanique Appliquee, vol. 4, University of Brussels, 1957, 
pp. 155-167. 

15. Mack, L.M.: Boundary Layer Linear Stability Theory. Special 

Course on Stability and Transition of Laminar Flow. AGARD-R-709, 
AGARD, Paris, France, 1984, pp. 3-1 and 3-81. (Avail. NTIS, 
AD-A147243 . ) 

16. Mack, L.M.: Review of Linear Compressible Stability Theory. On 

Stability of Time Dependent and Spatially Varying Flows, 

D.L. Dwoyer and M.Y. Hussaini, eds.. Springer, 1987, pp. 164-187. 

17. Moore, F.K.: Unsteady Laminar Boundary Layer Flow. NACA 

TN-2471 , 1951. 


31 



18. Moore, F.K.: Aerodynamic Effects of Boundary-Layer 

Unsteadiness. Sixth Anglo-American Aeronautical Conference, 

J. Bradbrooke and V.A. Libby, eds.. Royal Aeronautical Society, 
Folkestone, 1957, pp. 439-476. 

19. Patel, M.H.: On Laminar Boundary Layers in Oscillatory Flow. 

Proc. Roy. Soc. London Series A, vol. 347, no. 1648, Dec. 9, 
1975, pp. 99-124. 

20. Pedley, T.J.: Two-Dimensional Boundary Layers in a Free Stream 

Which Oscillates Without Reversing. J. Fluid. Mech., vol. 55, 
pt. 2, Sept. 26, 1972, pp. 359-384. 

21. Rosenhead, L. , ed.: Laminar Boundary Layers. Oxford Univ. 

Press, 1963. 

22. Ryzhov, O.S.; and Zhuk, V. I.: Internal Waves in the Boundary 

Layer with the Self-Induced Pressure. J. Mec., vol. 19, no. 3, 
1980, pp. 561-580. 

23. Stewartson, K.: On the Impulsive Motion of a Flat Plate in a 

Viscous Fluid. Q. J. Mech. Appl. Math., vol. 4, June 1951, 
pp. 182-198. 

24. Stewartson, K.: The Theory of Laminar Boundary Layers in 

Compressible Fluids. Oxford University Press, 1964. 


32 








.5 


FI6URE 5. - REAL AND IMAGINARY DISTRIBUTIONS OF g FOR M^ V 2; 
A = (1 + i )/2; AND c = 1 - 1/M^. 



J I I I I I 

i 1.0 1.5 2.0 2.5 3.0 


a 



(a) (b) 

FIGURE 6. - REAL AND IMAGINARY DISTRIBUTIONS OF g FOR M^ = >/2; 
A = (1 + i)/M; c = 1 - 1/M_. 
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